Vortex-line liquid phases: longitudinal superconductivity in the lattice London model 



T.J. Hagenaars 1 , E.H. Brandt 2 , R.E. Hetzel 1 '*, W. Hanke 1 , M. Leghissa 3 -**, and G. Saemann-Ischenko 3 
1 Institutfiir Theoretische Physik, Universitdt Wiirzburg 
Am Hubland, 97074 Wiirzburg, Germany 
2 Max-Planck-Institut fiir Metallforschung, Institut fur Physik 
Postfach 800665, 70506 Stuttgart, Germany 
3 Physikalisches Institut III, Universitdt Erlang en- Number g 
Erwin- Rommel- Strasse 1, 91058 Erlangen, Germany 

We study the vortex-line lattice and liquid phases of a clean type-II superconductor by means of 
Monte Carlo simulations of the lattice London model. Motivated by a recent controversy regarding 
the presence, within this model, of a vortex-liquid regime with longitudinal superconducting coher- 
ence over long length scales, we directly compare two different ways to calculate the longitudinal 
coherence. For an isotropic superconductor, we interpret our results in terms of a temperature 
regime within the liquid phase in which longitudinal superconducting coherence extends over length 
scales larger than the system thickness studied. We note that this regime disappears in the mod- 
erately anisotropic case due to a proliferation, close to the flux-line lattice melting temperature, of 
vortex loops between the layers. 



I. INTRODUCTION 

The effect of strong thermal fluctuations on the flux 
lines is an important aspect of the physics of high-T c 
superconductors in the mixed state. It is believed that 
these fluctuations melt the ground-state flux-line lattice 
into a vortex-line liquid below T c |l],|2| . 

In Nelson's analogy of vortex lines with world lines of 
2D bosons || , the absence of superfluidity of the bosons 
corresponds to the presence of superconducting coher- 
ence along the field direction in the vortex-line system 
|^,^). The liquid phase without longitudinal supercon- 
ductivity (sometimes referred to as an "entangled vortex 
liquid") corresponds to the 2D boson superfluid state. 
Similarly, the vortex-line liquid with longitudinal super- 
conductivity ("disentangled vortex liquid") corresponds 
to the normal liquid state of the 2D bosons. The quan- 
tity corresponding to the inverse 2D-boson temperature 
in the vortex-line system is the system size L z along the 
field direction. Therefore, for a small enough L z , Nelson's 
analogy predicts longitudinal superconducting coherence 
in the vortex liquid. In Rcf. J5|] it has been argued that 
the disentangled vortex liquid exists as a true thermo- 
dynamic phase (L z — > oo), at least for the case where 
the distance between the vortex lines is (much) smaller 
than the magnetic penetration depth \\. The scenario 
in which the longitudinal superconductivity is destroyed 
already within the vortex-lattice phase, leading to a "su- 
persolid" phase, has been considered in ||]. 

There exists a substantial recent literature that deals 
with MC simulations of the vortex- line phases, using vari- 
ous model approximations P,[7|-^6[ . In two recent papers, 
in which the vortex-line system was studied within lattice 
London models P,pT|, opposite conclusions have been 
reached regarding the presence of longitudinal supercon- 



ductivity within the vortex-liquid phase for systems with 
L z ~ 15 — 30 lattice constants. Although in both works 
cited, the longitudinal coherence was argued to exist only 
over a finite length scale (and therefore no true disentan- 
gled vortex- liquid phase was claimed to exist), a differ- 
ence of two orders of magnitude in the results for this 
length scale translates into opposite predictions for the 
experimentally observable behavior in typical samples 
studied in recent flux-transformer experiments, for tem- 
peratures above but close to melting. One of the motiva- 
tions of our work has been to clarify the reason(s) for the 
above disagreement. On the one hand, Chen and Teitel 
, who perform the Monte Carlo simulation at constant 
magnetic induction B, found that longitudinal supercon- 
ductivity persists far into the vortex liquid for L z = 30 
lattice constants. These authors argue that the length 
scale beyond which the vortex lines become entangled, 
corresponds to approximately 410 lattice constants. For 
YBCO this length scale for longitudinal coherence was 
estimated in Q to be much larger than the thickness 
of the s ample s studied in recent flux-transformer exper- 
iments [g!||30|]. As a result, their findings suggest the 
existence of a temperature regime in these experiments, 
in which the vortex lattice has melted but in which the 
flux motion is still maximally ^-correlated, in agreement 
with the results obtained for twinned YBCO samples [^9) . 
On the other hand Carneiro |l4| ] found evidence for the 
disappearance of the vortex-liquid regime with disentan- 
gled lines already between L z — 6 and L z — 12. These 
results were argued to agree with later experiments on 
untwinncd YBCO samples, in which the flux-line lattice 
melting and a loss of maximally z-correlated vortex mo- 
tion was found to coincide. 

In Carneiro's calculations the longitudinal response is 
calculated in a different way, discussed in more detail 
below, in which fluctuations in the net vorticities (in 
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two directions) are allowed. In a recent Comment p7[ , 
Chen and Teitel have criticized Carneiro's calculation^ 
scheme. They suggested that it effectively measures the 
transverse instead of the longitudinal response. Below, 
we will present results of MC simulations of the lattice 
London model. We present a direct comparison of the 
two different calculational schemes, and will point out 
that the method of Ref. |l4| has some serious problems. 

From our results we conclude that the isotropic lat- 
tice London model predicts a temperature regime in 
which longitudinal coherence over long length scales ex- 
ists within the vortex-liquid phase. We also present re- 
sults for a moderately anisotropic system, in which the 
temperature at which longitudinal coherence is lost, is 
roughly equal to the melting temperature. 

The outline of this paper is as follows. In section II 
we introduce the model and the Monte Carlo method(s). 
In section III we present our results. We summarize our 
findings in section IV. 



II. LATTICE LONDON MODEL AND MC 
ALGORITHMS 

The first MC study of the lattice London model was 
performed by Carneiro, Cavalcanti and Gartner jl£[ ]. The 
model was explored further by Carneiro for a supercon- 
ductor with free surfaces [|l9| . The lattice London model 
is formulated directly in terms of the vortex degrees of 
freedom. The vortex lines are modeled as consisting of 
elements of unit length based at dual lattice sites of a 
cubic lattice. At every dual lattice site we define three 
integers q^iR) — 0, ±1,±2, the vorticities in the di- 
rections fi = x,y, z. When a vorticity q^(Ri) is nonzero, 
its magnitude gives the number of flux quanta carried 
by the associated vortex-line element. The q fl (Ri) are 
subject to the continuity constraint 
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that ensures that the vortex-line elements form either 
closed loops or lines that end at the boundaries. Here 
Ri — e M runs over nearest neighbor sites of Ri. At con- 
stant B, the Hamiltonian of the isotropic lattice London 
model, expressed in terms of the vorticities, is given by: 



H = 4tt 2 J J2 q^Ri)q fl (Rj)9 l ,(Ri - Rj 



(2) 



where J = <&§<i / '(32ir 3 A 2 ) and g^{R) is the London inter- 
action with Fourier components 
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where k 2 = 2 — 2cosfc Al (k^ — 27m A1 /L M , 
0, 1, ...,L M — 1) and k 2 = J2p, K fi- Here we assumed an 
L z lattice with periodic boundary conditions. 
R = A 2 /A| = 1/r 2 is the anisotropy parameter, with 
Ai and A3 the magnetic penetration depths along the 
x-y plane and the z axis, respectively, and d the lat- 
tice constant. The Hamiltonian (H) can be either derived 



from the discrete version of the London free energy JL8 19 
or by a duality transformation of the Hamiltonian of a 
Villain-type lattice superconductor p0| . 

Monte Carlo sampling of the phase space for the vari- 
ables q^R) at constant B is performed as follows. The 
initial configuration is prepared to contain the number 
of vortex lines we want to study, depending on the value 
chosen for the magnetic induction B. A Monte Carlo up- 
date step consists of adding at a given site a closed d x d 
square loop of unit vorticity with an orientation chosen 
randomly from the six possible ones. The addition of a 
closed loop preserves the constraint (|J). The standard 
Metropolis algorithm is employed to accept or reject the 
new configuration. Obviously, when only closed loops are 
added, the magnetic induction B with components 



d 2 V 



(4) 



(V = L x L y L z ) will be constant throughout the simula- 
tion. 

To probe superconducting coherence one can consider 
the helicity moduli introduced by Chen and Teitel 
T M = lim,t v ^o T M (fc l/ ). Here T ^ L (k v ) is the linear response 
coefficient between a perturbation 8A ext (k v )e^ of the ex- 
ternal vector potential and the induced supercurrent 
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where k ■ A ext = (London gauge). For our lattice Lon- 
don model these moduli are given by 
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where (11, v, a) is a cyclic permutation of (x,y, z), 



q<j{kv) — qa{kv,kfj 



0), and (..) denotes a 



thermal average. We will focus on the quantity j z = 
lim^^o "f z (k x ). We note that for a uniaxial supercon- 
ductor with B parallel to the z axis, one has j z = 
limfc y ^o lz(ky) from symmetry, where r y z {k y ) is defined 
by the expression (^) with the indices x and y inter- 
changed. If 7 2 = 1, one has a perfect screening of the 
corresponding perturbation of the vector potential, and 
thus longitudinal superconducting coherence. A sharp 
jump in 7 Z from unity is interpreted as loss of longitu- 
dinal coherence. In Ref. |2|] it has been shown that j z 
will be non-zero as long as the vortex-line system retains 
a finite shear modulus at finite wave vector. Physically, 
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this may be understood by the finding of Ref. |22] that 
helical instabilities of the vortex lines driven by a longi- 
tudinal current, can be stabilized either by pinning or by 
a finite shear modulus. 

In the MC calculations, 7 M is estimated from the quan- 
tity 7 z (fc„) at small k v . In Ref. |j) this was done by ex- 
trapolation of the results for {q a {k v )q a {—k v )) to k v = 
using a fit to a polynomial in n 2 . It turns out that such 
a fit does not work for higher values of A and/or T than 
considered in Q . In Ref. (l5) this has been demonstrated 
for A = 12, r = 1. However, this problem disappears if 
one examines the full ratio that defines 1 — 7 M (£v) (see 
Eq. rtq )). This quantity has only a weak dependence on 
fc„ |23~| , so that one already obtains a reliable estimate for 
1 — j v from the data point for the smallest nonzero wave 
vector, k v = 27r/L„. This direct way of estimating the 
response was employed before in Ref. Jl^|. We illustrate 
it with an example in Figs. 10 and 11. In particular, 
1 — r y z (k x ) is almost independent of k x (when B is along 
z). Therefore, it is possible to obtain meaningful results 
for j z already for systems with L x equal to only a few 
times the average distance between the vortex lines. 

In Ref. |0, a different set up was introduced for 
measuring longitudinal superconducting coherence. This 
original set up starts from the observation that at a free 
boundary parallel to the y-z plane, superconducting co- 
herence in the z direction corresponds to Meissner shield- 
ing of an infinitesimal uniform applied field in the y di- 
rection 22 1 . For a lattice London model in the form of 
a slab, with fbc in the x direction and pbc in the other 
directions, and in a constant applied magnetic field, this 
shielding is measured by the transverse magnetic perme- 
ability (Jtry. 



dH y lH y=° 
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where (..} is now a thermal average weighted by the 
Hamiltonian 
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H = H-?—<Po/d 2 Y,UK 
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Here / M = H^d 2 /$□ is the dimensionless applied mag- 
netic field in units of the elementary flux quantum $o 
per plaquette area. In (||) f ^ multiplies the microscopic 
magnetic field 6 M : 



d 2 V 



(9) 



When fiy is zero, the external field H y is shielded by su- 
percurrents running parallel to e z and thus the system is 
a longitudinal superconductor. When fi y is nonzero, this 
is interpreted as a loss of longitudinal superconductivity. 
Note that = 1 — r y z (k x = 0), which is identically zero 
in the case of periodic boundaries and constant B. For a 
slab placed in a constant applied magnetic field however, 



the vorticity fluctuations measured by the thermal aver- 
age in the right hand side of Eq. (Q) can be nonzero by 
virtue of the free boundaries. At such boundaries the el- 
ementary vorticity fluctuations satisfying the constraint 
(0) are elementary loops with one side missing (the side 
outside the sample), so-called incomplete loops. 

The MC method at constant applied field, introduced 
by Carneiro in Ref. fll4| , uses the vortex- vortex interac- 
tion g(R) for a system with periodic boundary condi- 
tions in all directions, but updates at sites in the planes 
x = 1 and x = L x as if these planes were (adjacent to) 
free boundaries, allowing the addition of closed elemen- 
tary loops parallel to the y-z plane and of incomplete 
loops parallel to the x-y and x-z planes. The periodic 
boundary conditions are imposed also in the x direction 
(through the form of g{R)) to diminish finite size effects. 
Although x = 1 and x — L x are therefore no genuine 
free boundaries, Carneiro has claimed that the addition 
of incomplete loops incorporates the essence needed for 
the calculation of fi y . It is however not a priori evident 
whether or not this treatment of the boundaries leads to 
meaningful results for \i y . In fact, we think that it does 
not, as we will argue from a consistency check presented 
in paragraph III. A. 

In Carneiro's method, the allowed incomplete loops do 
not lead to net vorticity fluctuations in the x direction. 
We note that this allows us to simultaneously extract -f z 
from a simulation in which we calculate /i y , namely from 
the fluctuations of q x at k x = k z = and the smallest 
nonzero k y . As a result, we are able to compare the re- 
sults of the two methods in one and the same simulation. 



III. RESULTS 



We have considered vortex-line systems corresponding 
to B = 1/8, B = 1/15 and B = 1/24 (B along z and 
measured in units of $>o/d 2 ). The corresponding ground- 
state vortex lattices are shown in Fig. 1. 



(a) 



(b) 



(c) 



FIG. 1. Ground state vortex lattices, projected onto the 
x-y plane, for different values of the dimensionless magnetic 
field /: (a) B = 1/8, (b) B = 1/15, (c) B = 1/24. 

For each temperature we start in the ground-state con- 
figuration. Runs consist of 16384 MC sweeps through the 
lattice. Of these the first 8192 sweeps are discarded for 
equilibration, and the second 8192 are used to compute 
the thermal averages. 
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A. Isotropic case. 

We first present results obtained employing Carneiro's 
method, for different values of the magnetic field f x = 
f y = 0J Z = f. For B = 1/8 and B = 1/24 we chose 
the commensurate lattice size L x = L y — L z = 12 and 
for B = 1/15 the sizes L x = L y = 15, L z = 6,15. In 
Fig. 2 we show the vortex-lattice correlations for vortex 
elements q z {Ri) in the z-direction: 



P(R) = 



q z (Rj + R)q z {Rj)) 



(10) 



as a function of temperature, for B = 1/8, Ai = 12; they 
are in quantitative agreement with the corresponding re- 
sults reported in Refs. |l4| and jl5|]. We measure the 
temperature in units of the coupling J/ks- The melting 
of the flux- line lattice is reflected in the decay of the corre- 
sponding correlations. From the curves in Fig. 2 we may 
estimate the melting temperature to be T m = 2.7 ±0.1. 
This estimate is actually an upper bound of the true melt- 
ing temperature. As has been noted in Refs. P|Jl6|l, the 
discrete mesh introduces an artificial pinning potential 
for the vortex lines that tends to increase the apparent 
melting temperature. This effect is reduced when consid- 
ering smaller vortex densities, that correspond to a finer 
effective mesh. In ]T(| the depinning from the mesh was 
argued to occur below the (true) melting temperature for 
B lower than approximately B c = 1/32. 




•— • P(0,4,0) 
B — a P( 1,2,0) 
♦ — ♦ P(2,2,0) 
^ P(2,4,0) 
< — * P(2,6,0) 
▼ — ▼ P(4,0,0) 
P(4,3,0) 
P(0,0,0) 



incomplete elementary loops that were generated in the 
planes x — 1 and x = L x . 
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FIG. 3. Snapshot of the vortex- line configuration in the 
liquid phase for B — 1/15, Ai = 15, L x — L y — L z — 15, 
T = 2.4. The perspective is from above, the upper and lower 
x-y planes are enclosed by thin lines, and the fat lines indicate 
the vortex lines. Note that periodic boundary conditions are 
imposed in all directions, and that at some places different 
parts of one and the same vortex line connect through these 
boundaries. Due to the introduction of incomplete elemen- 
tary loop fluctuations in the planes x — land x = 15, the 
microscopic field is allowed to fluctuate (see text). 




T/T 



FIG. 2. Vortex-vortex spatial correlations P(R) for B — 
1/8, Ai = 12, L x = L y = L z = 12 as a function of tem- 
perature. Temperature is measured in units of J/fcs. The 
filled symbols represent correlation functions that measure the 
translational order corresponding to the ground-state config- 
uration depicted in Fig. 1(a). 

As an illustration of the liquid phase, in Fig. 3 a 
snapshot of the vortex configuration for the B = 1/15, 
Ai = 15, and T = 2.4 = 1.2T m is shown. In this snap- 
shot we observe a few (vortex-line segments built up of) 



FIG. 4. The transverse response \i y as a function of tem- 
perature, for the parameters (a) B = 1/15, Ai = 5, L m = 
L y = 15, L z = 6, (b) B = 1/8, Ai = 12, L x = L y = L z = 12, 
(c) B = 1/15, Ai = 15, L = 15, (d) B = 1/15, Ai = 5, 
Lx = L y = L z = 15, and (e) B = 1/24, Ai = 12, L = 12. 
Sample error bars are shown. For all curves the temperature 
is normalized to the melting temperature T m . These melting 
temperatures are T = 1.35 ± 0.1 (a) and T = 2.7 ± 0.1 (b), 
T = 2.15 ± 0.1 (c), T = 2.0 ± 0.1 (d), T = 1.7 ± 0.1 (e), 
and have been estimated from the decay of the vortex-lattice 
correlations as shown in Figs. 2 and 3 [25]. 
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In Fig. 4, curve (b) shows fj,y ELS Bi function of tem- 
perature for B = 1/8, Ai = 12. As in Refs. @|^, fx y 
is non-zero above the melting temperature. Curve (d) 
corresponds to B — 1/15, A = 5, the case studied at con- 
stant B in Ref. M. Here \i y is still small at the melting 
temperature, but rises fast for slightly higher tempera- 
tures. For L z — 6 (curve (a)) we find that \x y stays 
close to zero in a substantial temperature regime above 
T m . As has been noted before in Ref. Jig ], this shows 
that n v does not correspond to the transverse response 
1 — j x , as was suggested in Ref. |l7]] , because the latter 
response is resistive in the liquid phase. For further com- 
parison, we have included results for B = 1/15, A = 15 
(c) and B = 1/24, A = 12 (e). For all cases (b)-(e), with 
L z = 12 or L z = 15, fi y starts to rise at around the 
melting temperature. This is in sharp contrast with the 
results reported for y z in Refs. |@. 

We have also extracted j z from these same simulations, 
in the way explained in Sec II. In Fig. 5 we compare 
fly (triangles) and 1 — j z (circles) for B = 1/15, Ai = 
5, L z = 15, the same parameters that were considered 
in Ref. Up to T — 2.5T m , these results for j z are 
within error bars identical to those for "f z calculated us- 
ing the algorithm of Ref. |J] (shown as squares). We find 
that 1 — 7z rises from close to zero to close to one at a 
temperature approximately twice as large as T m . Simi- 
larly, for all other curves in Fig. 4 with exception of curve 
(b) (/ = 1/8), we find that j z is close to one up to the 
highest temperature values shown. Therefore, the tem- 
perature at which j z drops from unity lies well above the 
melting temperature, in agreement with the results ob- 



tained in Refs. [§ and Q. For B = 1/8 (curve (b) in 
Fig. 4, A = 12) the temperature region between the drop 
of P(R) and the drop of 7 Z is smaller. This is probably 
caused by the pinning artefact due to the discrete mesh, 
that tends to increase the apparent melting temperature 
more for B = 1/8 than for the lower vortex densities 
studied here 01. 




FIG. 5. Comparison of the responses fi y (triangles) and 
1 — 7z (circles) for B = 1/15, Ai = 5, L x — L y — L z = 15 as a 
function of temperature. Temperature is measured in units of 
J/fcs. The data points indicated with squares are results for 
1 — 7z from an independent simulation using the MC method 
with constant B. Sample error bars for 1 — 7 Z refer to the 
error in the data point for the smallest nonzero k v (see text). 




FIG. 6. Responses fj, y , 1 — y z , specific heat Ch, and absolute vorticities for B = 0, Ai = 12, L x — L y = 24, L z = 12. Sample 
error bars are shown for fi y , 1 — j z , and Ch- The insets show snapshots of the vorticity configurations at T = 4.0 and T = 5.7. 
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FIG. 7. Vortex-vortex spatial correlations P(R) for R — 
0.04, B = 1/15, Ai = 5, L x = L y = 30, L z = 10 as a function 
of temperature, together with 1 — j z . 




FIG. 8. Vortex-vortex spatial correlations P(R) for R — 
0.04, B = 1/15, Ai = 5, L x = L y — L z — 15 as a function 
of temperature, together with fj, v , 1 — 7 Z , and the absolute 
vorticities in x and y direction. Sample error bars are shown 
for fiy and 1 — <y z . 
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FIG. 9. Snapshot of the vortex- line configuration in the 
liquid phase for R = 0.04, B = 1/15, Ai = 5, L x = L y = 
L z = 15, T = 0.25. 



Summarizing the above results, we find that fi y gen- 
erally tends to rise at significantly lower temperatures 
than 1 — -f z . We will now present a consistency check 
on these response coefficients. We do this by consider- 
ing the B = superconductor-to-normal transition that 
is driven by (the unbinding of) thermally excited vortex 
loops. In this transition, the loss of superconducting co- 
herence is accompanied by a peak in the specific heat and 
a sharp increase in the number of vortex loops pq | , and 
the results for the response coefficients should therefore 
match these features. 

In Fig. 6 we present the results for fi y following 
Carneiro together with the helicity modulus ^ z follow- 
ing Chen and Teitel, calculated in the same simulation 
for an isotropic system at B = 0. We find that \i y rises 
at much lower temperatures than 1 — j z . The peak in the 
specific heat, an independent signature of the transition 
from the superconducting to the normal state, is situated 
around T ~ 5.9 and coincides with the rise in 1— j z . This 
provides strong evidence that j z and not fi y is the correct 
quantity to measure the superconducting coherence. In 
Fig. 6 we also show the absolute vorticity per unit cell as 
a function of temperature, which is a measure of the num- 
ber of thermally excited vortex loops. The increase in the 
vorticity is strongest close to T « 5.9, which is another 
signature of the transition p7|Jl5|| . The rise in fi y takes 
place at temperatures where the number of vortex-loop 
fluctuations is still very small and where superconductiv- 
ity can therefore not yet be destroyed. As insets in Fig. 6 
we show snapshots of the vortex loop configurations for 
T = 4 and T = 5.7. We conclude that fluctuations in the 
net vorticity q y that are induced by the allowed incom- 
plete loops in the planes x = 1 and x = L x and that lead 
to a nonzero /i y , do not signal a loss of superconducting 
coherence in the bulk |2Sj |. 



B. Anisotropic case. 

We now discuss results for a moderately anisotropic 
case in a non-zero applied field. In Figs. 7 and 8 we 
compare results for R = 1/25 (or T = 5), correspond- 
ing to YBCO. Here we took B = 1/15 and Ai = 5. 
The curve with data points shown as stars in Fig. 7 is 
the helicity modulus "f z following Chen and Teitel, or 
more precisely j z (k x ) at the lowest k x . The results for 
l—j z (k x ) are shown in Fig. 11. The lattice size was taken 
to be L x = L y = 30, L z = 10. We observe that for this 
anisotropic case, 1— j z rises to unity approximately at the 
vortex-lattice melting temperature. Similar results were 
obtained in a recent work by Nguyen et al. |]lq| , where 
the effect of anisotropy was studied systematically. Due 
to the anisotropy, the energy cost of vortex-loop fluctua- 
tions parallel to the x-y planes is reduced. As explained 
in Ref. p6fl , these parallel loops, which are most easily 
nucleated near a vortex line running in the z direction, 
contribute to the decoupling of the planes, so that lon- 
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gitudinal coherence is lost at lower temperature than for 
the isotropic case. In Ref. jl6| it was shown in addition 
that for strong anisotropy R < 1/100, the parallel loops 
proliferate already at a temperature (far) below the melt- 
ing temperature of the vortex lattice, thus giving rise to 
a vortex lattice phase without longitudinal coherence in 
between these temperatures. 




FIG. 10. Transverse response 1 — ~/ x (k y ) for different tem- 
peratures, calculated at constant B, for R = 0.04, B — 1/15, 
A = 5, L x = L y = 30, L z = 10. The peak for higher T is at a 
wave vector commensurate with the vortex lattice. 




FIG. 11. Transverse response 1 — -y z (k x ) for different tem- 
peratures, calculated at constant B, for R = 0.04, B = 1/15, 



A = 5, L x = L y = 30, L z = 10. 



In Fig. 8 we show Carneiro's permeability \i y calculated 
at the same parameters as considered in Fig. 7. Here we 
used L x = L y = L z = 15. We observe that there is a 
similar reduction of the temperature at which longitu- 
dinal coherence is lost when compared to the isotropic 
case. However, \i y rises already well below the melting 
temperature. The only physical mechanism that could 
be responsible for this early rise would be a proliferation 
of parallel loops. However, there is no such prolifera- 
tion at this temperature, as can be seen from the snap- 



shot of the vorticity configuration at T = 0.25 shown 
in Fig. 9. There are only few parallel loops present, as 
can be seen from the absolute vorticities shown in Fig. 8 
as well. However, according to the result for fj, y , the 
planes should be almost decoupled at this temperature. 
This again gives strong evidence that \i v is governed by 
artefacts. In Fig. 8 we also show 1 — 7 Z calculated si- 
multaneously from the fluctuations of q x at finite k y . In 
agreement with the calculations shown in Fig. 7, this 
longitudinal response is found to turn resistive at the 
vortex-lattice melting temperature. 



IV. SUMMARY 

We investigated the occurrence of longitudinal super- 
conductivity in the vortex liquid phase in a clean type-II 
superconductor, within the lattice London model, mo- 
tivated by recent contradicting results reported in Refs. 
0j and Jli) . For an isotropic superconductor, our results 
show three distinct regimes as a function of temperature. 
Between the Abrikosov lattice phase and a vortex liquid 
without longitudinal coherence, there is a liquid regime 
in which longitudinal superconducting coherence extends 
over length scales larger than the system thickness stud- 
ied Q. In the context of flux-transformer experiments, 
this intermediate regime should correspond to resistive 
regime with z-correlated vortex motion for samples that 
are sufficiently thin. In the moderately anisotropic case 
R = 1/25 (corresponding to YBCO), we find that lon- 
gitudinal superconducting coherence is lost close to the 
flux-line lattice melting temperature. This is due to a 
proliferation, close to the melting temperature, of vortex 
loops between the layers |l6| ]. In recent flux-transformer 
work |30| on YBCO single crystals, flux-line lattice melt- 
ing and a loss of maximally ^-correlated flux motion was 
found to coincide. 
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